What do we learn from the shape of the dynamical susceptibility of glass-formers? 
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We compute analytically and numerically the four-point correlation function that characterizes 
non-trivial cooperative dynamics in glassy systems within several models of glasses: elasto-plastic 
deformations, mode-coupling theory (MCT), collectively rearranging regions (CRR), diffusing de- 
fects and kinetically constrained models (KCM). Some features of the four-point susceptibility Xii^) 
are expected to be universal: at short times we expect a power-law increase in time as t* due to 
ballistic motion (f^ if the dynamics is Brownian) followed by an elastic regime (most relevant deep in 
the glass phase) characterized by a t or \/f growth, depending on whether phonons are propagative 
or diffusive. We find both in the /3, and the early a regime that ~ ^''i where is directly related 
to the mechanism responsible for relaxation. This regime ends when a maximum of X4 is reached at 
a time t — t* oi the order of the relaxation time of the system. This maximum is followed by a fast 
decay to zero at large times. The height of the maximum also follows a power-law, X^i^*) ~ ^*^- 
The value of the exponents fj. and A allows one to distinguish between diff'erent mechanisms. For 
example, freely diffusing defects in rf = 3 lead to /.i = 2 and A = 1, whereas the CRR scenario 
rather predicts either /i = 1 or a logarithmic behaviour depending on the nature of the nucleation 
events, and a logarithmic behaviour of X4(^*)- MCT leads to /i = b and A = I/7, where b and 7 are 
the standard MCT exponents. We compare our theoretical results with numerical simulations on 
a Lennard- Jones and a soft-sphere system. Within the limited time-scales accessible to numerical 
simulations, we find that the exponent jj, is rather small, fi < 1, with a value in reasonable agreement 
with the MCT predictions, but not with the prediction of simple diffusive defect models, KCMs with 
non-cooperative defects and CRR. Experimental and numerical determination of Xii't) for longer 
time scales and lower temperatures would yield highly valuable information on the glass formation 
mechanism. 

I. INTRODUCTION 

The idea that the sharp slowing down of super-cooled liquids is related to the growth of a cooperative length 
scale dates back at least to Adam and Gibbs m. But it is only a few years back that this idea has started being 
substantiated by convin cing exp eriments 1 3 , B , I M JHl , numerical simulations 0, S S 0j 01 simple 

microscopic models E El Ea El E ETIzl ll^, HI HI IH - One of the basic problem has been to find an 
observable that allows one to define and measure objectively such a cooperative length scale. An interesting quantity, 
proposed a few years ago in the context of mean field p-spin glasses |27j (see |28j for an important early insight) and 
measured in simulations, is a four-point density correlator, defined as 

G4(r,t) = {p{0,Q)p{0,t)p{r,Q)p{r,t)) - (p(0, 0)p(0, t)) (p(f, 0)p(f, t)), (1) 

where p(r, t) rep resents the density fluctuations at position r and time t. In practice one has to introduce an overlap 
function w |27j . to avoid sin gula rity due to the evaluation of the density at the same point or consider slightly 
different correlation functions |29j. This quantity measures the correlation in space of local time correlation functions. 
Intuitively if at point an event has occurred that leads to a decorrelation of the local density over the time scale t, 
■ G4{r,t) measures the probability that a similar event has occurred a distance r away within the same time interval 
■ ■ ' t (see e.g. [lOl)- Therefore G4(r, t) is a candidate to measure heterogeneity and cooperativity of the dynamics. The 
best theoretical justification for studying this quantity is to realize that the order parameter for the glass transition 
is already a two-body object, namely the density-density correlation function C(t) = (p(0, 0)p(0, t)), which decays to 
zero in the liquid phase and to a constant value in the frozen phase. The four-point correlation Gi{f, t) therefore plays 
the same role as the standard two-point correlation function for a one-body order parameter in usual phase transitions. 
Correspondingly, the associated susceptibility X4(t) is defined as the volume integral of Gi{r,t), and is equal to the 
variance of the correlation function (27. 31, 33|. The susceptibility Xii'^) h^-s been computed numerically for different 
model glass formers, and indeed exhibits a maximum for t = t* ~ t^, the relaxation time of the system pdl Il2 . Il3l [l4l | . 
The peak value X4:{t*) is seen to increase as the temperature decreases, indicative that the range of G4(r, t*) increases 
as the system becomes more sluggish. The dynamical correlation length ^4(t*) extracted from Gi{f,t*) in molecular 
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dynamics simulations grows and becomes of the order of roughly 10 intcr-particlc distances when the time-scales is 
of the order of 10^ microscopic time-scales tq with tq ~ 0.1 ps for an atomic liquid. In experiments close to the glass 
transition the dynamical correlation length has been found to be only slightly larger, between 10 and 20 inter-particle 
distances 0, 0|. This is puzzling because experiments are done on systems with relaxation times that are several 
orders of magnitude larger than in simulations. In fact, extrapolating simulation results in the experimental reg ime 
would lead to much larger dynamical correlation lengths. The origin of this puzzle is still unclear, see Ref. |l8l | for 
a recent discussion. Experiments on dynamical heterogeneity bridging the gap between numerical and macroscopic 
time-scales would be extremely valuable to resolve this paradox. 

Several scenarii have been proposed to understand the existence of non-trivial dynamical correlations, and their 
relation with thermodynamical singularities. Adam and Gibbs 0, Kirkpatrick, Thirumalai, Wolynes and collabo- 
rators [S^l (for a different formulation, see Ref. H^l), and Kivelson and Tarjus [S^ have proposed, using somewhat 
different arguments, the idea of collectively rearranging regions (CRR), of size ^ that increases as the temperature 
is decreased. The evolution of the system is such that these regions are either frozen or allowed to temporarily and 
collectively unjam for a short time until a new jammed configuration is found. 

In apparent contradiction with the existence of the growing length scale, the Mode Coupling Theory (MCT) of 
glasses states that the self consistent freezing of particles in their cages is a purely local process with no diverging 
length scale at the transition |36l |. However, this point of view is in disagreement with the results found for mean- field 
disordered systems jl^H^ that are conjectured to provide a mean-field description of the glass transition and display 
an MCT- like dynamical transition. Indeed it was recently shown that within MCT G4{f,t) in fact develops long 
range correlations close to the critical MCT temperature Tc [sJl . Within a phase-space interpretation of the MCT 
transition, the mechanism for this cooperative behaviour for T > is the progressive rarefaction of energy lowering 
directions [s^l . Within a real-space interpretation, the MCT transition is due to the formation of a large number of 
mctastable states, each one characterized by a surface tension that increases from zero at Tc. As one approaches Tc 
from above, the relevant eigenvectors of the dynamical Hessian become more and more extended, which means that 
the modes of motion that allow the system to decorrelate are made of very well defined, collective rearrangements of 
larger and larger clusters of particles (see the recent work of Montanari and Semerjian [Sg). For smaller temperatures, 
T < Tc, 'activated events' are expected to play a crucial role. They are believed to be responsible for the destruction 
of the freezing transition at Tc- This regime has been tentatively described by adding 'hopping terms' in the MCT 
equations [3^ or within a CRR scenario |U . 

Exploiting yet a different set of ideas, models of dynamical facilitation, such as the Frederickson- Andersen 
or Kob- Andersen models Q j have recently been proposed as paradigms for glassy dynamics 0, l^^, . In these 
models, the motion of particles is triggered by 'mobility defects' that diffuse and possibly interact within the system. 
As the temperature is lowered or the density is increased, the concentration of defects goes down, and the relaxation 
time of the system increases. Dynamics is obviously heterogeneous since it is catalyzed by defects that cannot be 
everywhere simultaneously. The characteristic length scale in this case is related to the average distance between 
defects to some model and dimension dependent exponent 0, 0, |23, |U . 

Understanding the mechanism behind the growth of the dynamical correlation length is certainly an important 
step — arguably the most important one — to understand the cause of the slowing down of the dynamics. Furthermore, 
the different scenarii for the glass transition can be tested contrasting their quantitative prediction for the four-point 
correlation function G4(f, t) to the numerical, and hopefully soon experimental, results. Following these premises we 
investigate in this paper the analytical shape of G4(r, i) for several simple models. We show that G4{f,t) indeed 
contains some important information concerning the basic relaxation mechanisms. However we show that, perhaps 
disappointingly, models where cooperativity is absent or trivial lead to four-point correlation functions and dynamical 
susceptibilities Xi that exhibit non trivial features. Other, more complex observables will have to be defined to really 
grasp the nature of the collective motions involved in the relaxation process of glasses |^ [s^ . 

Let us summarize the main results of our study, in terms of the susceptibility X4(^) and time-sectors. In a super- 
cooled liquid there are separate regimes of time-scales corresponding to different physical behaviour (sec Fig. 1). 
On microscopic time-scales particles move ballistically if the dynamics is Newtonian, or diffusively if the dynamics 
is Brownian. On a longer time-scale, interactions start playing a role, which can be described approximately using 
elasticity theory, before a truly collective phenomenon sets in. This non trivial glassy regime is the /^-regime, within 
which correlation functions, as for example the dynamical structure factor, develop a plateau. The /3-regime is 
divided further in an early and a late /3-regime corresponding, respectively, to the approach and the departure from 
the plateau of the correlation function. Finally the structural relaxation time-scale on which correlation functions 
decay to zero is the a-regime. All previous studies have focused on the behaviour of X4 (0 times of the order of 
which correspond to the peak of X4(^)- We show that Xi{t) has in fact a rich structure in time and different behaviour 
in different time-sectors. In many of these regimes, x^i^) behaves as a power-law of time, f^, with different values of 
fi. During the ballistic time-scale one finds = 4 (/i = 2 for Brownian dynamics) whereas during the elastic regime 
(most relevant deep in the glass phase), the exponent becomes /-t = 1 for ballistic phonons and /x = 1/2 for diffusive 
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FIG. 1: Sketch of the time behaviour of X4(^)i with all the different time regimes, within the MCT description that we find to 
be a good description around Tc. As the temperature is lowered, we expect the elastic regime to extend up to r^. 

phonons. The behaviour in the /3 and a regimes is intimately related to the physical mechanism for relaxation and 
indeed we find quite different answers depending on which scenario we focus on. MCT predicts exponents /x = a and 
fi = b on time-scales corresponding respectively to the early and late (3 regimes, where a and b are the standard MCT 
exponents obtained from the study of the dynamical structure factor. The power-law t'' extends until the peak in 
X4{t) is reached. 

The other scenarii only make predictions in the a regime. In the case of CRR one has Xi ^ t or X4 ^ (Int)'^^^^'^' 
before the peak depending whether one assumes that the relaxation occurs via bulk nucleation events or domain 
wall fluctuations, see below. For diffusing defects in dimension = 3, the exponent is /i = 2. If defects have a non 
trivial diffusion exponent z, such that their displacement at time t scales as t^^^, then /i = 2d/z for d < z and ii — 2 
otherwise. The overall behaviour of Xi{t) is summarized by Fig. ^ which specializes to the MCT predictions for 
simplicity. 

Another important feature of X4 is the growth of the peak compared to the growth of the time t = t* ^ Ta at 
which the peak takes places [23|. This is found to scale as X4(t*) t*"^, with A = (logarithm) for CRR, A = 1 for 
freely diffusing defects, A = d/z for anomalously diffusing defects for d < z and A = 1 again for d > z. Note that if 
the defect diffusion coefficient itself scales with t* as , as for example in the one-spin facilitated FA model, there 
is an extra contribution that gives \ = 1 — fiord>z. Finally, one has A = I/7 in the context of MCT, where 7 
describes the power-law divergence of the relaxation time as the critical MCT temperature is approached. 

We have checked these predictions in two model systems of glass- forming liquids: a Lennard- Jones and a soft-sphere 
mixture. Concerning the behaviour of Xii't) ii^ the late P and a regime, the most interesting time-sectors, we have 
found reasonable agreement with the MCT predictions for four point correlators. This agreement is by no means 
trivial and is actually quite unexpected unless MCT indeed captures some of the physics of the problem. Instead 
models of diffusing defects do not describe well the numerical results unless one assumes anomalously diffusing defects 
with z substantially smaller than 2. This is perhaps not very surprising since we are focusing on two fragile liquids 
(at least in the numerical time window) at temperatures well above the experimental glass transition. It might be 
that the predictions of these models work only on larger time-scales. In any case, we expect instead that for strong 
liquids displaying an Arrhcnius behaviour the predictions for Xii^) obtained studying model of simple diffusing defects 
should hold quantitatively, since it is indeed quite well established that relaxation in strong liquids is triggered by 
the diffusion of connectivity defects ^3 ■ Finally, the CRR picture does not agree quantitatively with our present 
numerical data. However, this picture is supposed to describe the liquid dynamics precisely in the low temperature/ 
long time regime that is presently beyond numerical capabilities. Again, experimental results probing the behaviour 
of Xi{t) in this regime would be highly valuable to put strong constraints on the different theoretical scenarii of glass 
formation. 

The organization of the paper is as follows. In Section^we discuss the behaviour of X4(t) on microscopic time-scales. 
Then, we analyze the predictions of elasticity theory in Section UTTI In Sections ll VI and Ivl we focus on the behaviour of 
X4{t) in the (3 and a regimes for MCT and CRR. In Section iyil we discuss the predictions of defect models analytically 
using an independent defect approximation and by numerical simulations of kinetically constrained models. In Section 
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IVIII we compare the different predictions to the results of numerical simulations of models of glass-forming liquids. 
We present our conclusions in Section IV 1111 



II. MICROSCOPIC DYNAMICS 



On very short time-scales the behaviour of Xi can be computed exactly. For simplicity, we characterize the dynamics 
through the self-intermediate scattering function, 



(2) 



and define the dynamic susceptibility as the variance of the fluctuations of Fs{k, t): 

2v 



l^cosfc-irKO-rKO)] 




(3) 



The full intermediate four point scattering function defined in E g. flll in fact contains very similar information, even 
for interacting systems - as shown by numerical simulations 0, Il3j|. 

On a very short time-scale particles move ballistically if the dynamics is Newtonian, fi{t) — fi{Q) = Vit + 0{t^), 
where Vi is the velocity of the particle i at time t. Since the system is in equilibrium all the Vi^s are independent 
Gaussian variables with variance {vi ■ Vj) = SijSkBT/m, where T is the temperature, m the mass of the particles, and 
ks the Boltzmann constant. Using this property it is straightforward to obtain 



Fs{k, t) = exp —k 



knT . 



2m 



and 



Xi{t) = Fs{k,tf 



cosh ^ 
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For an interacting particle systems this is only valid on short time scales, for example smaller than the collision time 
for short ranged interactions. This leads to an initial power-law increase that reads 
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(6) 



Note that if one had chosen Langevin dynamics (i.e. dtfi — dpH + fji) instead of Newtonian dynamics, Eqs. H5I6(I 
would have been identical except for the replacement of ksTt^ /m by 2Tt, again for small times. Thus changing from 
Newtonian to Langevin dynamics, the initial power-law increase of X4(t) changes from to . This is similar to 
the change in the mean square displacement that increases as and t, respectively, for Newtonian and Langevin 
dynamics. 

In the above example, however, it is clear that the increase of X4 with time has nothing to do with the increase of a 
correlation length, since particles arc assumed to be independent. In other words, the four-point correlation G4{r,t) 
has a trivial (5-function spatial dependence, but the height of the 5 peak increases with time. As will be discussed 
later in the paper, it is important to normalize Xii^) by the value of Gi^f = 0,t) to conclude from the four-point 
susceptibility that a length scale is indeed growing in the system. 



III. ELASTIC CONTRIBUTION 



For longer time-scales the interaction between particles starts playing a role. Gcnerically one expects that in the 
time regime where the displacements of particles remain small, an elastic description should be valid. In a solid, or 
in a glass deep below Tg, there is no further relaxation channels and the elastic contribution to X4 should be the only 
relevant one. In a super-cooled liquid around the Mode-Coupling temperature Tc, the elastic regime is interrupted 
by the collective P regime, where in some sense phonon-phonon interactions completely change the physical picture. 
Although we expect such a crossover, we have at present no detailed theoretical description of it. 
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In the following we analyze again the behaviour of the four point self-intermediate scattering function assuming 
that the dynamical behaviour of the liquid can be described, within a restricted time-sector, as an elastic network 
(we will discuss later how to include, in a phenomenological way, viscous flow). Perhaps surprisingly, we find a non 
trivial structure for G^i in this model, with an ever growing 'cooperative' length scale which comes from the dynamics 
of phonons, that represent the simplest form of cooperativity. 

We consider an isotropic solid immersed in a viscous thermal bath. The energy of the system is given by: 

where ki, K2 are the Lame coefficient, Uij — \[^^ + is the deformation tensor and (p the displacement field from 

an undeformed reference state. Note that is simply the continuum limit of the displacement of each particle with 
respect to its equilibrium (bottom of the well) position. 

As is well known, the above energy leads to three independent phonon modes (one longitudinal and two transverse 
modes). For simplicity, we only consider one deformation mode and write the Hamiltonian in Fourier space as: 

where k is an effective elasticity modulus. The mode k has an energy Ek = nk^ (f)k(t>~k/'^ and therefore we expect, 
in equilibrium, ((/)k4'-k) = T/nk"^, where the Boltzmann constant has been set to unity. Our goal is to calculate the 
dynamical correlation functions of the system. We describe the dynamics by a Langevin equation with a local noise: 

m^^ + .^ = «A0(f,t)-fC(.^,O, (9) 

where C,{x, t) is a Gaussian noise uncorrelated in space and time, of variance equal to 2uT. Taking the Fourier 
transform: 

"^-^ + i^-^ = "Kfc^^fe + Cfc(i) (10) 

Cfc(t) is again a Gaussian noise uncorrelated for different fc's and time. 

In this section, we only consider in details the over-damped case m = and set D = n/v, but also give at the 
end the result for the purely propagative case v = (see also Appendix A). One easily deduces the non equal time 
correlation in the over-damped case: 

(0.(t)^_,(O)) = ^e-^^^*. (11) 

Let us now define the function: 

F(')(r,t) = ^<5(r-r,(0))cos(9h(t)-r,(0)]), (12) 

i 

whose average equals the self- intermediate scattering function up to a constant (the particle density) . 
Using the microscopic definition of (p we obtain that 

C{q,t) = {F^''\r,t)) ~ {e'i[HP,t)-^im]) = g-g^([(0(r,t)-0(r,o)]=>/2^ ^^3) 

where the last equality comes from the Gaussian nature of the deformation field. Using the above results on the 
correlation of the Fourier modes, we find: 

(W.V)-*KO)P> = f /l^^ij (14, 

As is well known, this integral behaves differently for d < 2 and for d > 2, reflecting the fact that phonons destroy 
translational order in low dimensions. As above, we will only consider here the physical case d — 3, relegating the 
discussion of the other cases to Appendix A. For d = 3, we need to introduce an ultraviolet cutoff A on the wavevector 
k, which is the inverse of the underlying lattice spacing a. Then, the above integral goes to a constant oc A at large 
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times, reflecting the fact that particles are localized in their 'cage'. Therefore, the self-intermediate scattering function 
C{q,t) decays at small times h?Dt <C 1 before saturating to a 'plateau' value given by: 

fq = C{q, oo) = exp f-c^-^) , (15) 



where c is a numerical constant. [Note that TAq^ /k has no dimension, and is expected, from a Lindemann criterion, 
to be of the order of 0.05 at half the melting temperature and for q = A]. In real glass- forming liquids, this plateau 
phase does not persist for ever, and C{q,t) finally decays to zero beyond t = Tq, in the so-called a-relaxation regime. 
A modification of the model to account for this decorrelation will be discussed later. Furthermore, the above pseudo-/? 
regime predicted by elastic theory does not explain quantitatively the /?-regime in super-cooled fragile liquids, except 
probably on relatively short time scales, say up to a few picoseconds. On the other hand, at temperatures below Tg 
or for strong glasses, we expect that the elastic regime will extend up to Tq and compete with other mechanisms, such 
as the defect mediated correlation discussed in Section IVll below. 

The calculation of G^^\f,t) = {F^'^'> {r' ,t)F^'^'> {r' + r,t))c is detailed in Appendix A. One immediately sees that 

G^\f,t) is governed by a diffusive correlation length ^(t) ^ ^/Dt with D = k/v, as expected from the structure of 
the Langevin equation that describes relaxational dynamics. Clearly, in the case of propagative phonons, one finds 
^ Vt with = K/m. The final result, see Appendix A, is: 

Gi'\r,t) = C\q,t) {cosH2q^R{r,t)) - l) , (16) 

where 

R{r,t)^-{Dt)^-''/'-F(^] (17) 



K \VDt, 

and we find (see Appendix A) F{z) ~ (inz)^^ for z ^ 1 and F{z) ~ (27r^/^)~^ exp(— z^/4)/z^ for z ^ 1. Note 
the similarity between the expression in (|16|l and the corresponding one jSJ derived in the previous section. One can 
check that indeed the short time behaviour is indeed the one derived before in the case of Langevin dynamics for 
the particles, as expected. Let us now focus on long-times, but still within the elastic regime: A^Dt ^ 1, and for 

G^\f,t)^fl[cosH^^)-iy (18) 

Suppose for simplicity that we are in a regime where the argument of the cosh is always small, corresponding to 
the limit Tq^A <C k (remember that by definition r > a = 2tt/A, where a is the inter-atomic distance). Then, 
G4{f,t) ^ for <C r <C £,{t). For larger scales r ^ ^{t) decays as a Gaussian, i.e. super-exponentially fast. 
Note that the small r behaviour of Gi{r^t) is not of the Ornstein-Zcrnike form (l/r in d = 3). Integrating G4 over r 
we find the dynamical susceptibility, 

xi'\t) ^ —i^at)- (19) 

This result is actually valid both for in the diffusive limit where £^{t) = \/Dt and in the propagative regime where 
i{t) = Vt. Therefore XiHt) increases either as or as t (note that in the limit of small times one recovers the or 
laws obtained in the previous section). In the general case, one expects a crossover between a propagative regime 
at small times t < m/v = D /V"^ (of the order of ps in glass formers, see ^iid a diffusive regime for longer time 
scales. Thus, looking at x'f' (0 ^ function of time in a log-log plot one should see first a straight line corresponding 
to the ballistic or diffusive motion leading respectively to slope /i = 4 or = 2, bending over towards a smaller slope 
(1 or 1/2, or both depending on the strength of the dissipation). The order of mag nitude of xi'^^(t), as given by 
Eq. H19|l . can be estimated to be ^ 10~^ — 10~^a^^(t) for q~ A. In the propagative regime with < = 1 ps, y = 3 10^ 

m/s, a = 0.3 nm, one finds ^ = 10a and x^i^ ~ 10"^ - 10"^ , i.e. a small, but perhaps detectable signal from the 
phonons. Only on much larger time scales will the elastic contribution be significant, a regime that can be reached 
deep in the glass phase 01 ■ mentioned above, other collective modes come into play in super-cooled fragile liquids, 
in particular around the mode-coupling temperature, and give rise to the /3-regime where 'cages' themselves become 
more complex, extended objects [SJ. 

The above calculation shows that in an elastic solid with diffusive or propagative phonon modes, the dynamical 
susceptibility increases without bound, rcfiecting the presence of Goldstonc soft-modes in the system. Of course, in 
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a real glass the correlation function C^'^^{t) eventually decays to zero beyond the a relaxation time Tq, as particles 
start diffusing out of their cages, far away from their initial position. If phonons were the only relevant excitations, 
this would cause the dynamical susceptibility to peak around t ^ t* ^ t^. A phenomenological model that describes 
the decay of x'f' (t) within the above elastic framework is to assume a (Maxwell) viscoelastic local modulus: 



dt 



dt' 



(20) 



with 7 ^ T~^, corresponding to a frequency dependent elastic modulus k{lu) = iKUj/{iuj + 7). In this model, the 
dynamics of 4> becomes diffusive at times > 7~^, and the dynamic structure factor therefore decays exponentially 
beyond that time. Of course, the model itself becomes inconsistent at large times, since the underlying lattice needed 
to define the deformation field </> has by then totally melted. 

The conclusion of this section, however, is that since super-cooled liquids behave at high frequencies (uj ^ 7,t^^) 
like solids, the four-point correlation and dynamical susceptibility are expected to reveal, in a certain time domain, a 
non trivial behaviour unrelated to the structure of the 'collective processes' discussed below (MCT, diffusive defects, 
CRR's) that one usually envisions to explain glassy dynamics. 



IV. MODE COUPLING THEORY 



As mentioned in the introduction the mode-coupling theory of supercooled liquids predicts the growth of a coopera- 
tive length as the temperature is decreased or the density increased (27.„ .28.. .3L| , and makes detailed predictions on the 
shape of Xii^)- The four-point correlation function becomes critical near the mode-coupling transition temperature 
Tc- The behaviour of the susceptibility X4(0 encoded in ladder diagrams [28II3H . From the analytical and numerical 
results of [sj, and analyzing the ladder diagrams we have found that in the /3 regime: 

Xiit)-^ fi{te'^'n/Ve t^T0 (21) 

and in the a regime 

X4(i) ^ /2(ie^/'"+^/''')/e tr^T^ (22) 

where e = T — Tc, a, b and 7 = l/2a -I- 1/2& are the MCT exponents for the dynamical structure factor, and /i(x) 
and f2ix) are two scaling functions. Requiring that the dependence on e drops out when te^/^° <C 1 one finds that 
/i(.t) ^ x°' when a; <C 1. This leads to a power-law behaviour, X4 ^ in the early (3 regime, i.e. when the 
intermediate scattering functions approaches a plateau. In the same way, matching the behaviour of /i when a; ^ 1 
to the one of /2 when x 1 one finds another power-law behaviour, ^4 ~ t*", on timescales between the departure 
from the plateau and the peak of X4- We give in Fig. ^a schematic summary of the shape of Xiit) within the MCT 
description of supercooled liquids. 

Finally, as discussed in |3ll |. at times t ^ t* ^ Ta, Xi reaches a maximum of height (T — Tc)~^. Using the relation 
Ta - (T - Tc)-'', valid within MCT, one finally finds x^t*) - T^/'^. 



V. COLLECTIVELY REARRANGING REGIONS 



Under the term CRR, we gather many similar scenarii that differ in their details, as discussed in the introduction 
P, lU I33. 135I I . Within the frustration-limited domains scenario of Ref. [s^ it seems natural to envision the dynamics 
as the activated motion of domains pinned by self-generated disorder. In the case of the random first-order theory of 
Refs. 113 13 1 the details of the decorrelation mechanism are not entirely clear. There should be, on the one hand, 
activated fluctuations of domain walls between different states, again pinned by self-generated disorder. However, 
the fluctuations leading to a change of state may be the nucleation of a completely different state starting from the 
bulk. The latter process can be modeled as a nearly instantaneous event with a certain (small) nucleation rate. In 
the following we shall analyze separately these two types of fluctuations and their consequences on the shape of Xiit)- 



A. Instantaneous events 



Suppose that the dynamics is made of nearly instantaneous events that decorrelate the system in a compact 'blob' 
of radius ^o- The probability per unit time and volume for such an event to appear around site r is F. We compute 
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the four-body correlation of the persistence, nr{t), defined to be equal to one if no event happened at r between times 
and and equal to zero otherwise. The four-body correlation is then defined as 

G4(f, t) = {nrit)no{t)) - (23) 

Clearly, the averaged correlation function, C{t) = {nr{t)), is simply given by C{t) = exp(— fiF^Qi) where is the 
volume of the unit sphere. For G4{f,t) to be non zero, an event must have happened simultaneously at r and at 0, 
leading to 

G4(f, t) = C^t) [exp (Fteo'/(r/eo)) - l] , (24) 

where f{x) is the volume of the intersection between two spheres of unit radius with centers at distance x apart. 
Clearly, f{x > 2) = 0. Therefore, Gi{r,t) is non zero only if r < 2^0, and is in fact roughly constant there. 
For a given r satisfying this bound, first grows linearly with time, reaches a maximum for t = t* ^ F^^^q' 
and decays exponentially beyond that time. The same behaviour is found for Xi{i)^ that grows initially as with 
/i = 1, and reaches a maximum such that X4(<*) oc ^q. Assuming finally that these events are activated Is^j . 
with a barrier growing like T^q, where tp is a certain exponent, one expects t* ^ tq exp(T^Q'/T), and therefore 
XA{t*) cx (Int*)'^/'^ (X Co- 

The rearranging regions could have of course more complicated shapes than the simple sphere assumed above. As 
long as these objects arc reasonably compact, the above results will still hold qualitatively. On the other hand, if these 
regions are fractal with a dimension df < d/2, the above results on G^ will hold with the argument in the exponential 

given by Ttr'^'^f^'^] one also finds t* « I/F^q^ and xtit*) ■ 

B. Domain wall fluctuations 

In this case the picture that we have in mind is similar to the case of a disordered ferromagnet with pinned domain 
walls, where the typical time to flip a domain is comparable to the inter-event time. In that case, an 'event' is in fact 
the slow fluctuation of domain walls that progressively invade the bulk of the domain. The early time behaviour of 
t) is given by the square of the number of particles that relax per unit volume thanks to the same domain wall (see 
for the same situation out of equilibrium in pure systems). Let again be the typical size of a domain and £(t) 
the lengthscale over which the domain walls fluctuate during time t. Considering that on the surface of each domain 
there are order {£,o/£)'^~^ subdomains of linear size £ and that the number of particles in each of these subdomains is 
proportional to t"^, we get Xti^) oc ^q'^{^o/()'^^^£'^'^ oc £'^'^^/^o- We are descarding for simphcity both the possibihty 
of fractal domains and that transverse fluctuations behave differently from longitudinal ones. Assuming thermal 
activation over pinning energy barriers that grow like T£'^ we finally get Xiit) oc S,Q^{lntY~^^^^ . Therefore, in 
this case, the exponent fj, is formally zero and the growth of X4(^) is only logarithmic. The maximum of X4 occurs at 
time t* such that £{t*) w ^o, which implies that the maximum of the susceptibility also scales logarithmically with t* , 
X4{t*) OC £,Q^{lnt*)'^^^^^ oc Cq. The same scalin g of the maximum of the susceptibility with the typical domain size is 
obtained in non-disordered coarsening systems [30j . 

The conclusion of the above analysis is that if the CRR relaxation is due to both instantaneous events and domain 
wall fluctuations, the latter will dominate the time behaviour of Xi before the peak as can be readily deduced by 
comparing their relative contributions to Xii't)- If for some reason, domain walls are particularly strongly pinned and 
bulk nucleation becomes dominant, then the exponent /i = 1 should be observable. The height of the peak, on the 
other hand, behaves identically in both models. Thus, as the temperature is reduced, one should see a power-law 
behaviour before the peak with an exponent < 6 < 1 in the MCT regime followed by an effective exponent /i either 
decreasing towards zero or increasing towards one depending on whether the domain wall contribution dominates or 
not. However, at lower temperatures, the elastic contribution will also start playing a role, that might completely 
dominate over the CRR contribution. This suggests that other observables, that quantify more specifically the 
collective dynamics, should be devised to reveal a CRR dynamics. 

VI. DEFECT MEDIATED MOBILITY 

A. Independently diff'using defects 

As the simplest realisation of the defect mediated scenario for glassy dynamics advocated in [Tsl IT^ ITqI l20l , 
we consider a lattice model in which mobility defects, or vacancies, perform independent symmetric random walks. 
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We assume for the moment that these vacancies cannot be created or destroyed spontaneously. We shall compute the 
same function 6*4 (f,t) as in Eq. (|23|l above, arguing that when such a vacancy crosses site r, the local configuration 
is reshuffled and the local correlation drops to zero. Therefore, nr{t) is equal to one if no vacancy ever visited site r 
between t = and <, and zero otherwise. Thus, {nr(t)) represents a density-density dynamical correlation function 
whereas {no{t)np{t)) — {no{t))^ corresponds to Gi{r,t). 

From now on we will denote by the number of vacancies, by V the total volume, by pv = Ny/V = 1 — p the 
vacancy density and by P§{t) the probability that a vacancy starts in z at time zero and never reaches x till time t. 
The probability that a vacancy starts in z at time zero and reaches for the first time x at a, time u < t is therefore 

The computation of {nr^{t)) is identical to the target annihilation problem considered in l46i|. Since we assume 
defects to be independent, the defect distribution is uniform and we have: 



cxp 



-Pv - Pv 



E 



p. it) 



(25) 



The correlation function {nx{t)ny{t)) can be also expressed in terms of probability distributions of a single random 
walk in a similar way: 



{n^{t)ny{t)) 



^ E 



1 - 



z,z^x,y 

2 



1 ^ {l^P^^(t)~Pl^{t) 



z^z^x^y 



V V 



V E ^^^W-^ E iPy,^it) + Pl-yit)) 



z,z^x,y 



z,z^x,y 



exp -2py - p, 



Z.Z^X 



p.pm 



Pv_ 

2 



E iPydt)+P^^^{t)) 



(26) 



z^z^x.y 



where P^-(t) is the probability that a vacancy starts in z at time zero and never reaches neither x nor y till time t, 
P^-.{t) is the probability that a vacancy starts in z at time zero and reaches x sX u < t but never reaches y till time t. 
In Eqs. H25I26|I we are left with the calculation of probabilities of the form P^it), P^^yit) + Pyxit) for a single random 
walk. This can be done using Laplace transforms and, concerning P^{t), the computation has been performed a while 
ago the details can be found in Appendix B. 

In the continuum limit, {x — y)/ \J Dt/2 ~ 0(1), i.e. for independent Brownian motion with diffusion coefficient Z?, 
the final expression for {nx{t)) on time scales much larger than one is, in three dimensions, 



(ri^(t)) = exp[-/9,„ - ciDpyt] 



(27) 



where ci is a constant fixed by the short-lengthscale physics, i.e. the underlying lattice structure (see Appendix B). It 
is clear from this expression which is valid in all dimension larger than two that the relaxation time-scale is governed 
by the vacancy density p„ and reads r = 1/ {cipyD). Physically r corresponds to the time such that each site has 
typically been visited once by a defect. 

The final expression for G4 is, for time and length scales much larger than one, and in the small vacancy density 
limit, — > 0, 



Pv 



■ exp 



2t 



2 nU g~ 2D t 

du I dv- -—7- 

Jo {2Dvtf/^ 



(28) 



where C2 is a constant of order unity. Note that the correlation length at fixed t is given by ^(i) = \/ Dt. For r <C (,{t), 
G4(r, t) ^ 1/r whereas for r 3> ■^(i), G4 decays at leading order as a Gaussian, that is, much faster than exponentially. 
The 1/r behaviour is cut-off on short-length scales, where H28|) does not hold. For r = one finds, when t ^ 1, 



Gi{r = 0,t) = {nx{t)) - {nx{t)f = exp(-t/T) [1 - exp(-t/T)] , 



(29) 



which behaves as t/r at small times. 
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By integrating (|28|l over r we get the dynamical susceptibility, 

X4W = ^(^)'cxp(-f). (30) 

For short times, t < t, the dynamical susceptibility is proportional to t^, so that fj, = 2. This is due to the diffusing 
nature of the defects. The main contribution to X4 is given by the square of the number of sites visited by the same 
defect, which behaves as pv{Dt)^ = -j- (|r) , since a random walk in three dimensions typically visits t different sites. 
For i > r, on the other hand, the correlation decreases because sites start being visited by different vacancies. The 
maximum of X4(<) is reached for t = t* = t, for which one has Xi{t*) ^ Pv^ ^ Dt* . Note that because random 
walks are fractals of dimension df ~ 2, the above relation can also be written as X4(i*) ^ a'^^'^f £^^f {t*), where we 
have added the lattice spacing a to give to Xi the dimension of a volume. If for some reason D depends on py, as it 
happens for example for the one-spin facilitated FA model where D (x py, then one finds t* ^ Py^ and Xiit*) t*^/^. 
Taking the Fourier transform of G'4(r, t) given by Eq. H28|l . we find the four point structure factor, S^^k, t), 

S4{k, t) = Xi{t)T{Dk^t)- T{u) = ^{u-l + e-") . (31) 

Note that S{k = Q,t) = Xii't)^ it should. Furthermore for large and small k, S4{k,t) behaves respectively as 
S4 ~ k~^ and 6*4 ~ X4 + O(k^), just as the Ornstein-Zernike form, though the detailed k dependence is different. 

One can also study this problem in dimensions c? = 1 or d = 2. Qualitatively, the same conclusions hold (diffusive 
correlation length \/Di, correlation time t* set by the density of vacancies, etc.), although the quantitative results 
differ because a random walk in d < 2 visits a number of sites that grows sub- linearly with time, see Appendix B.l 
and B.3. One finds in particular that Xii^*) ^ {Dt*Y/^ ^ ^'^it*), with logarithmic corrections for d = 2. The above 
arguments can be generalized if for some reason the vacancies have an anomalous diffusion motion, in the sense that 
their typical excursion between time t = and time t scales as t^/^ , where z is the dynamical exponent. When z = 2, 
usual diffusion is observed, but many models like diffusion in random media or kinetically constrained models may 
lead to sub-diffusion, where z > 2 |23,|4g. In this case, one expects the small time behaviour of Xi{t) to be given by 
X4{t) ~ t^'^/^ for d < z and for d > z with logarithmic corrections for d — z. Similarly, the behaviour of Xii^*) is a 
power-law, Xii^*) with A = d/z for d < z and A = 1 for d > z. 

In the above model, mobility defects were assumed to be conserved in time. However, it is certainly more realistic 
to think that these defects can be spontaneously created and disappear with time. Suppose that defects are created 
with a rate F per unit time and unit volume, and disappear with a rate 7 per unit time. The equilibrium density of 
defects is then py = T /^. The above results on Xi can easily be generalized. At small times, the number of pairs of 
visited sites will now behave as py{Dt)'^ — ^r{Dt)^/D. Because of the death of vacancies there is an extra decay of 
the dynamical susceptibility. The dominant rate of decay depends on the adimensional number jt. 

A very similar model for glassy dynamics was suggested in [49l | , where free volume is described as a diffusing coarse 
grained density field p(r, t) with a random Langevin noise term. Mobility of particles is allowed whenever the density 
p exceeds a certain threshold pQ. The mobile regions are then delimited by the contour lines of a random field, which 
already gives rise to a quite complex problem of statistical geometry jSOf. The particle density correlation in this 
model is a simple exponential with relaxation time r ~ exp(po/p), where p is the average free volume density. One 
can also compute Xii^) in this model to find, in d = 3, 



exp( -^j h -cxpf-^ 



(32) 



which behaves very much like the point like vacancy model studied above, with in particular, X4(0 ^ for t <^ t. 

Let us finally note that from the point of view of interacting particles on a lattice we have studied the persistence 
dynamical susceptibility, instead of the density-density correlations discussed in the introduction. This is because for 
the lattice gas problem at hand, the former does not show any interesting properties: except when a defect passes by, 
the local state is always the same, i.e. occupied. For completeness, we give the corresponding results in Appendix 
B.4. In a real system, however, the local configuration is going to be affected by the passage of a mobility defect, 
and one can expect that the density-density correlations will in fact behave more like the persistence dynamical 
susceptibility computed before. The correspondence betwe en p ersistence and self-intermediate scattering function is 
studied explicitly in kinetically constrained models in Ref. 



en p e 



B. Kinetically Constrained Models: numerical results 



Kinetically constrained models (KCM) postulate that glassy dynamics can be modeled by forgetting about static 
interactions between particles, putting all the emphasis on dynamical aspects. Among those models are, for example. 
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the Frcdrickson- Andersen (FA) , or the Kob- Andersen (KA) on hyp er-cubic lattices [Tsl . l25l | . The dynamics of these 
models can be understood in terms of diffusion of defects 0, 121I l25j and the models can be classified into cooperative 
and non-cooperative models, depending on the properties of such defects. For cooperative models the size, the 
density, and the time-scale for motion of the defects depends on the particle density (for conservative models) or 
temperature (for non-conservative models) and change very rapidly with increasing density or decreasing temperature 
|25| . KA and FA models with more than one neighbouring vacancy needed in order to allow the motion of other 
vacancies belong to this class. On the other hand for the one-spin isotropically facilitated FA model, a single facilitating 
spin is a mobile defect at all values of temperature and the model is non-cooperative. A recent analysis pll | suggests 
that for these models defects can be considered as non-interacting in d > 4, while for d < 4 the role of fluctuations 
becomes important. Therefore we expect that the previous results for the independent diffusing defects model should 
apply exactly for FA one-spin facilitated in c? > 4. Furthermore, since the corrections to the Gaussian exponents are 
not very large |2l| in three dimensions, we still expect a semi-quantitative agreement. In particular the initial increase 
of the dynamic susceptibility as X4(t) ^ N{t)^, where N{t) is the total number of distinct visited sites, is expected to 
be quite a robust result. Also, we expect a diffusive growth of the dynamical length scale ^{t) governing the scaling of 
Gi, at least in the limit ^{t) ^ ^o- At smaller times, one expects a crossover between a CRR regime when Dt ^ 
(where the dynamics inside the defects becomes relevant in cooperative models to a mobility defect regime for longer 
times. Hence, in principle, looking at the detailed properties of 6*4 (r,t) one should be able to extract the defect 
properties: density, size, time-scale and decide which theoretical scenario is most consistent with numerical results. 

In the following, we discuss numerical results for the one-spin facilitated FA model both in d = 1 and d = 3, and 
for the d = 1 East model where facilitation is anisotropic ■ The two models can be described respectively in terms 
of diffusive and sub-diffusive non-cooperative defects and indeed the numerical results are in quantitative agreement 
with the predictions of previous section, as will be explained in detail. We do not address the case of cooperative 
KCM models, for which a more complicated behaviour is expected. Indeed a first slowing down of dynamics should 
occur near a dynamical crossover displaying the properties of an MCT-like avoided transition psj . In this regime 
the model cannot be approximated as a system of independent freely diffusing defects and deriving a quantitative 
prediction for the behaviour of four point correlation and susceptibility would deserve further work. Such avoided 
transition should then be followed at lower temperature or higher density by an asymptotic behaviour described in 
terms of cooperative diffusing defects. 



1. One dimension 



Let us start with the simplest model, the d = 1 FA model. For a given temperature, we consider the time evolution 
of the following quantities. The analog of the spatial four-point correlator for this model is 



^4 ('^. = ^ E [(^. (t)n,+r it)) - (t)] , (33) 

i=l 

where n{t) = -/V^^ X]i3=i("-i(*)) is the mean persistence, ni{t) being the persistence at site i. We also measure the 
corresponding four-point structure factor 



1 ^ 

S4{k,t) = ^ E [{ndtMt)) - n'it)] e*(^-'"), (34) 

and as usual we get the four-point susceptibility as the fc ^ limit of the structure factor, Xii^) — S4{k — 0,t). We 
generally find that the results are in good agreement with the free defect model described above, at least at sufficiently 
low temperatures. 

In Fig. 121 we show the evolution of the spatial correlator (|33|l at a given low temperature, T = 0.2, and various 
times. At this temperature, the relaxation time is about r ^ 10^, so that time-scales presented in Fig.|2cover a range 
of times both smaller and larger than r. The dynamic susceptibility Xii^) has the usual shape with a maximum at 
a time close to r indicating that dynamics is maximally heterogeneous there. This non-monotonic behaviour of Xi 
in fact does not show up in the spatial correlators of Fig. |21 which display instead a smooth monotonic evolution 
with time. The spatial decay of 6*4 (r,t) becomes slower when t increases indicating the presence of a monotonically 
growing dynamic lengthscale ^(t). 

One can estimate the time dependence of ^(t) by collapsing the data of Fig. |21 using a form like: 

G4(r,0 ^G4(0,t)g (^0 . (35) 
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r 

FIG. 2; Four-point spatial correlator 1331 in the d = 1 FA model at fixed temperature, T — 0.2, and various times, t = 10'', 
3.10^, 10", 3.10*, 10^ 10®, 3.10", 6.10® (from left to right). The correlator is normalized by its r = value. At this temperature, 
the relaxation time is r ~ 10®, so that time-scales cover both regimes where t/r is smaller and larger than 1. 

Doing so, we find that ^ ^ i"'^^ is a reasonable representation of the data at T = 0.2. Correspondingly, we find that 
the increase of Xii^) for t < t \s well-described by a power-law, Xi ~ t"'®^, so that the expected scaling Xi ^ 
is reasonably verified given the unavoidable freedom in estimating the range of time-scales where power-laws apply. 
The values of these exponents are not far from the ones expected from freely diffusing defects in one dimension, 
although slightly smaller. Indeed, we recall that the results in Appendix B.l. predict ^ = Vlji, Xii^) °^ pH^Y ^-iid 
Xi{t*) = Pi where p is the density of defects, D their diffusion coefRcient and t* the time at which Y4 (t) reaches its 
maximum value. This last prediction is also in good agreement with the numerical results (see e.g. pOl'). 

Repeating the simulation at lower temperature. T ~ 0.15, we obtain X4 ~ t^'^"^, showing that deviations from theo- 
retically expected values are partly due to preasymptotic effects that presumably disappear at very low temperatures. 

It is important to remark that the scaling form (|35() is only approximately supported by the data. The scaling in 
fact deteriorates when times become larger than t. This can be seen in Fig. |21 where data for large times become 
more and more stretched, indicating an increasing polydispersity of the dynamical clusters. Note that a change in the 
shape of the spatial correlator makes a quantitative determination of ^ problematic. Usually, one wants to collapse 
various curves using a form like H35|l to numerically extract ^. Strictly speaking, this is not possible here if one works 
at fixed T and varying t over a large time window. This difficulty provides a second possible explanation for the small 
discrepancy between the measured values of exponents and the theoretical expectations. 

The observation of a monotonically growing length begs the question: how can the correlation length increase 
monotonically with time while the volume integral of the spatial correlator x^ is non-monotonic, as reported in the 
previous section? This is due to the fact that we have presented in Fig. [21 results for the normalized correlator, 
Gi{r,t)/G4^{r = Q,t). By definition, 6*4(0, i) = n(l — n), hence the normalization itself exhibits a non-monotonic 
behaviour. If one considers the normalized susceptibility, X4 = [^4(6, t)N]^^ [(?^ff^m) — '^^(0] : o^^^ indeed finds 

that X4 is monotonically growing as well. 

In numerical works, the quantities that have been studied are in fact, most of the time, normalized, and the 
corresponding X4(^) observed for realistic systems shows a peak, at variance with what is observed in the d = 1 FA 
model. As we shall show below, this is due to the one-dimensional nature of the model, and this difference is not 
observed in three dimensions. This difference in the behaviour of the normalized dynamical susceptibility between 
one and three dimensions is indeed in full agreement with the independent defect diffusion computation, see previous 
section and Appendix B. 

Results are qualitatively similar in the one-dimensional East model. The dynamic susceptibility X4(0 develops a 
peak that grows and whose position is slaved to the increasing relaxation time when temperature decreases. At fixed 
temperature, a monotonically growing lengthscale is observed, while the scaling relation xa ~ still holds within our 
numerical precision. The novelty of this model lies in the fact that exponents arc now temperature dependent, as all 
other dynamic exponents in this model. For instance, we find that ^{t) ^ t°-^^ a.t T — 0.4, ^{t) ^ at T — 0.2. 
These results are in agreement with the above predictions of the independent defects model if the defect motion is 
sub-diffusive, with a dynamic exponent z = Tq/T, as expected from |l7j . Due to the quasi one-dimensional nature 
of the relaxation process in the three-dimensional generalization of the East model |l8j , these results most probably 
carry over to larger dimensions where they would differ by numerical factors only. 
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FIG. 3: Four-point correlations in the d — 3 one-spin facilitated FA model both in real space (left) and in Fourier space (right) 
at fixed temperature, T = 0.17, and various times indicated in the figures. In Fourier space, points represent numerical data, 
while full lines are fits to the form 1361 with fitting parameters described in the text. 

2. Three dimensions 

In d = 3, the situation is more subtle. Results for the normalized susceptibility of the one-spin facilitated FA model 
were presented in Ref. j22| |. where it was found to have the standard non-monotonic shape already described several 
times above. We find that the non-normalized x^i^) has the same qualitative behaviour. Therefore, contrary to the 
d = 1 case normalization is not a crucial issue in three dimensions. 

In the following we check the predictions for independent diffusing defects in three dimensions for the susceptibility 
and correlation length obtained above, i.e. ^{t) = Voi, Xii^) p S,{tY S'^d X4(^*) = where p is the density 
of defects, D their diffusion coefficient and t* the time at which Xii^) reaches its maximum value. We find a 
semiquantitative agreement with above prediction, with small deviations in the exponents that should be due to the 
interaction among defects. In particular the scaling of the peak with the density of defects was already analyzed in 
[2^ . where the result Xi(t*) ^ obtained, with e ~ 0.03. As for the correlation length, we find £,{t) oc t*^ "*^, 

that shows again a small deviation from the diffusive prediction. Regarding the increase at f ^ t of the susceptibility 
we find a power-law as predicted. As in d = 1, the exponent changes slightly when decreasing temperature because the 
scaling regime where power-law applies becomes more and more extended. We find Xi ~ ^^'^ at T = 0.25, Xi ~ i^'^^ 
at r = 0.17 and X4 ~ t^'^^ at T = 0.095. This seems to indicate that the deviation from the scaling Xiit) °c t'^ 
calculated for the independent diffusing defect model is partly due to preasymptotic efiiects that are less and less 
important at lower temperature. Unfortunately, we were not able to measure ^ at much lower temperatures with 
sufficient accuracy. We expect that even at very low temperature a small deviation from the exponent of independent 
defects should survive due to the interaction among defects. 

In Fig. Owe show the four-point correlations both in real and Fourier space, Eqs. H33|) and H34|) . In these curves the 
temperature is fixed at a low value, T = 0.17, and time is varied in a wide range that includes the relaxation time, 
t(T ~ 0.17) ~ 5.10^, where the dynamic susceptibility also peaks. For times t <^ t, the spatial decay of Gi{r,t) is 
fast. When t increases, the spatial decay becomes slower, once again indicative of an increasing dynamic correlation 
length £,{t). When t becomes larger than r, however, spatial correlations seem to become weaker. It is obvious from 
Fig. that the volume integral of Gi{r,t)/G/i{Q^t) decreases when t grows larger than r. This is very different from 
the one-dimensional case in Fig. |21 but consistent with all known numerical results. 

However, a closer look at Fig. |3| reveals that even though the initial spatial decay of 6*4 (r,i) is stronger at larger 
times, the contrary is true at large distances. This indicates that the topology of the dynamic clusters changes when 
t grows larger than r, but that ^(t) may keep increasing in a monotonic manner. Since the spatial correlator is very 
small at large distances, quantitative measurements of ^(i) are more easily performed in Fourier space via Si^k^t). 

At short time, a fit of S'4(fc, t) using the functional form given by Eq. H31|l works reasonably well, but the fit quickly 
deteriorates at long time. We have therefore used the following generalization of Eq. (|31|l : 

S4k,t)=Xiit)J^p[k^e{m J',3iu) = — {u-l + e-Y - (36) 

Freely diffusing defects correspond to /3 = 2 and £,{t) ~ ^/i. Using P{t) as an additional free parameter, we are able 
to fit S4{k,t) at all times, see Fig.|3| We find that (3 decreases from /3 w 2.5 at small times to /3 w 1 for the longest 
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FIG. 4: Time dependence of the dynamic susceptibility in the binary LJ mixture at two different temperatures. The lines are 
power-law fits with exponents indicated in the label. 



times scales investigated, which corresponds to t w 5t. At such large times, the dynamic susceptibility has already 
decreased by a factor w 300 from its maximum value at t = r, and correlations become very weak indeed. The values 
for (3 found from the fits are consistent with the value /? ~ 2.15 reported in Ref. |22| where only fixed time ratio 
t/T{T) = 1 at different temperatures have been studied. From this fitting procedure, we deduce a monotonically 
growing dynamic length £,{t), even beyond t = t{T). Fitting its time dependence with a power-law, we get ^ ^ t^''^'^ 
which appears to be slightly sub-diffusive, but close to the value found above in the one-dimensional case. 

In conclusion we find that on small enough time-scales, one indeed has good agreement with the above calculations 
based on freely diffusing defects, therefore defect branching and defect coagulation can be neglected. However, for 
longer time-scales, significant deviations appear which correspond to the evolution of the exponent /3(t) and should be 
responsible for the small deviations of the predicted exponent for Xi- Physically, the time evolution of the exponent 
(3{t) characterizing the large k behaviour of the dynamic structure factor is reasonable. At very short-times, dynamic 
clusters consist of coils created by random walk ers, and an exponent close to /3 = 2 can be expected. For times < ~ r, 
clusters look critical, as described in Refs. [22, |22|, and the exponent /3 = 2 — 77, 77 < is expected. At very large 
times, clusters are most probably extremely polydisperse because the remaining spatial correlations at large times are 
due to the largest regions of space that were devoid of defects at time and that take therefore a large time to relax. 
But at large times, some isolated sites that have not been visited by defects during the relaxation might survive so 
that the distribution of dynamic clusters at large times is very wide, see Ref. for snapshots. A small value of /? 
can therefore be expected. 



VII. NUMERICAL RESULTS ON ATOMISTIC MODEL SYSTEMS 



In this section, we study numerical results for the dynamic susceptibility and structure factor of a super-cooled 
liquid simulated by molecular dynamics simulations. The model we study is mainly the well-known binary Lennard- 
Joncs mixture as first defined and studied in Ref. [s^ ]. but we report also some results for a soft-spheres mixture 
studied in |3^.j5^j53 . We do not give details about our numerical procedures since these were given several times in 
the literature |2l[ 129115^ . 



A. Dynamical susceptibility 

In previous works on various realistic liquids, the dynamic susceptibility was reported several times [^ ll0llT^l27j . 
It is known to exhibit at peak at a time-scale slaved to the quantity chosen to quantify local dynamics. Typically, 
particle displacements are chosen, and one computes therefore the variance of some dynamical correlation, 

X4it) = N [{F^t)) - {Fit)r] , (37) 

with 

1 ^ 
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The dynamic quantity Fi{t) can be chosen as som e 'p ersistence' function in which case {F{t)) resembles the overlap 
function usually measured in spin systems [TolIT^ . l2^. Other choices are |2lll29| 

F,{t) = cosiq-Sn{t)), (39) 

where g is a wavector chosen in the first Brillouin zone, and Sfi{t) is the displacement of particle i in a time interval t. 
In the limit of small |fc|, it is better to study Fi{t) = \S'ri{t)\/y^ Ar'^{t), where Ar^(t) is the mean square displacement 
of the particles 0, 0| . 

Whereas the general shape of Xii't) is well-documented in the literature, its precise time dependence was never 
discussed. In Fig. 2) we present the time dependence of Xi (0 the binary Lennard- Jones mixture at two different 
temperatures. The data are presented in a log-log scale, in order to emphasize the existence of several time regimes 
that are generally hidden in the existing reports. To build these curves, we choose as the local observable, for a 
wavector that corresponds roughly to the typical inter-particle distance. 

In the ballistic regime at very short-times, we find that x^i^) ~ described from Sectional The system 

then enters the time regime where dynamic structure factors typically exhibit plateaus, as a result of particle caging. 
As seen in Fig 01 this is also the case for Xii^)- Finally, X4{t) reaches a maximum located close to the relaxation 
time extracted from the time dependence of {F{t)), and then rapidly decays to its long-time limit, equal to 1/2 
in the present case. In Fig. ^ we fitted the time dependence of the increase of Xii't) towards its maximum with 
power-laws Xi ^ The fits are satisfactory, although they only hold on restricted time windows. We find a slight 
temperature dependence of the exponent fi. For instance, we find fi w 0.9 at T = 0.47, and ^ w 0.73 at T = 0.42. As 
already discussed in the case of kinetically constrained models above, it is not clear how the restricted time window 
used to determine the exponents might afi'ect their values. However, the data in the Lennard-Jones system behave 
quantitatively very differently from both theoretical results obtained from freely diffusing defects and numerical results 
in the one-spin facilitated d = 3 FA model, where fi = 2. The small temperature evolution in the LJ liquid differs even 
qualitatively from the one-spin facilitated d = 3 FA model where the exponent was found to increase when decreasing 
temperature. These observations tend to discard a description of this super-cooled liquid via a scenario with simple 
independently diffusing defects, even interacting ones. The above value of fi is in principle compatible with the 
predictions of elasticity theory, which yields fi = 1/2 ot fi = I depending on the damping of phonons. However, the 
time scale in which the above mentioned power-law behaviour holds in the Lennard-Jones mixture corresponds to the 
/3 regime where the displacement of particles is no longer small and the elastic description unjustified. Within MCT, 
on the other hand, Xi should increase in that regime with an exponent fj, = b that is known from previous analysis of 
the dynamics of the binary Lennard-Jones mixture, b « 0.63 |52j |. The values found above are somewhat larger, but it 
is hard to know how preasymptotic effects influence the numerical data. Moreover, the value closest to 6, /i « 0.73, is 
obtained for T = 0.42, a temperature already lower than the mode-coupling singularity located at ~ 0.435 in this 
system. MCT also provides a prediction for the height of the peak, xt ^ t*'^/^ , where 7 was measured to be « 2.3, 
leading to A = I/7 ~ 0.43. This prediction is in good agreement with the results of Ref. |22| where Xiit*) ^ t*'^'^ was 
reported. 

If on the other hand one insists to use a non-cooperative kinetically constrained model to describe the Lennard- 
Jones liquid, the small value of the short time exponent /i forces one to choose a 'fragile' KCM model, such as the 
East model described above, where the exponent for the dynamic susceptibility is found to be much smaller than 
the diffusive value ^ — 2, and indeed to decrease when temperature is decreased. However, the explanation given 
in Ref. that the large dynamic length-scales observed in the Lennard-Jones system are due to the fact that the 
underlying KCM model is relatively strong becomes hard to reconcile with the present results. This findings can 
therefore be added to the list of unusual features displayed by supposedly 'fragile' numerical models for super-cooled 
liquids [S^. On the other hand, note that our results do not discard the possibility that cooperative KCM (in a 
proper density or temperature regime) display a four point correlation and susceptibility quantitatively similar to the 
one of the Lennard-Jones liquid. Indeed, as stressed in e.g for these models one expects a first regime of slowing 
down of dynamics due to an avoided mode-coupling transition. The susceptibility and four point correlation could 
then well be quantitatively comparable to that of Lennard-Jones liquids. 

Finally, it is of course a natural question to ask whether the above agreement between MCT predictions and 
numerical results is only restricted to the Lennard-Jones system. Using the unpublished data of Ref. 0| for a soft- 
sphere binary mixture where Tc ~ 0.22 — 0.24 113, 113 we actually found very similar results. Close to Tc a power-law 
behaviour of Xi a function of time can again be observed. For instance, Xi ~ t'^'^^ ^or T = 0.26. In Fig. [3 we 
plot X4, defined as in Ref. |23|, as a function of time. We also display the power-law behaviour predicted by MCT 
before the peak with the exponent b ~ 0.59 taken from Ref. |53l |. There is a similar agreement between the exponent 
A measured from the height of the peak and the value of I/7 extracted from an MCT analysis of the data. 

The fact that the predictions of MCT for the four-point susceptibility are in reasonable agreement with numerical 
simulations in both systems is significant, since the exponents b and I/7 are measured on (local) two point functions 
and /i and A on four-point functions. The relation between these exponents test a rather deep structural prediction 
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FIG. 5: Dynamic susceptibility y4(t ) at T — 0.3 and 0.26 (from left to right) in a log-log plot as a function of time for the 
soft-sphere binary mixture of Ref. I37ll54l . The data was kindly provided to us by D. Reichman and R. A. Denny. The straight 
line represent the MCT prediction for the power-law behaviour before the peak. 



of MCT that relates time scales to length seales [3l| . More numerical work, on other model systems with different 
values of 6, for example, would be needed to establish more firmly whether the coincidence observed in the present 
paper is or not accidental. 



B. A growing length scale? 



We focus now more directly on the dynamic lengthscale. In previous works, the dynamic lengthscalc ^ extracted 
from four-point correlations was measured either at fixed temperature for various times t where it was found to be 
non-monotonic [l^ ITsL but monotonic in or at fixed time t = t(T), for different temperatures, where it is 
found to be increasing when the temperature decreases 0,0,1^. In practice, to extract $,{t,T) from the four-point 
correlation function either in real space or in Fourier space, one needs to postulate a specific functional form of G4. 
In this respect, the results of the previous section on simple lattice KCM's with no underlying liquid structure prove 
instructive. It is clear that with data similar to Fig. |3 but obtained with much smaller system sizes, with much less 
statistics, and polluted by the underlying structure of the liquid, the precise extraction of dynamical length-scales from 
Molecular Dynamics simulations is not an easy task. More fundamentally, extracting ^ from fitting either G4(r, t) or 
S4{k,t) to a time-independent scaling form necessarily biases the data as discussed above. This also shows that it 
is a much easier and safer procedure to work say at i = t(T) and different temperature to observe the growth of a 
cooperative length ^(r, T) when decreasing T. On the other hand, it is not a priori granted that the growth law of ^ 
with t = t{T) when changing T is identical to that of ^(i, T) with t at a given temperature T. We will not be able to 
answer this question with our numerical data. 

With the above caveats in mind, we present in Fig. some numerical data in the binary Lennard-Jones mixture 
at a fixed temperature, T = 0.5, and three different times which fall before, at and after the peak in Xi{^)- The 
difficulty of getting clear-cut quantitative determinations for ^ are obvious from Fig. One would need much larger 
system sizes to properly measure S'4(fc, t) at small wavectors, large times and low temperatures. The system simulated 
here contains 1372 particles. One could possibly increase the number of particles by a factor 10, but the increase in 
linear size would be very modest, a factor 10^/'^ ~ 2.15. Nonetheless, we have fitted the data in Fig. with a simple 
empirical form, 

for < fc < fco, fco ~ 7.21 being the position of the first peak in the static structure factor. As for the d = 2> FA 
model, the exponent /3(t) and the dynamic length are fitting parameters. There is an additional free parameter, 
the additive constant G in Eq. H40|) . which accounts for the fact that the structure of the liquid starts to be visible 
and creates some signal in 54(fc, t) when k — s- fep. The results of the fitting procedure are presented in Fig.|Slwith lines 
going through the data. Note that the fits in Fig.^lare constrained at low k by the value of the dynamic susceptibility 
Xiit)- The most important result from Fig.|Slis that if the functional form of Si{k,t) is given some freedom, here via 
the time dependent exponent (3{t), the extracted dynamic lengthscale ^(i) indeed continues to grow monotonically 
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FIG. 6: Left: dynamic susceptibility at T = 0.5 and q = 4.21. The vertical lines indicate the times at which S4{k,t) is 
evaluated in the bottom figure. Right: the corresponding three S4{k,t) (the last two have been multiplied by 2 for clarity). 
Lines are fits to the form (I40II . the fc limit being fixed by the value of x-i{^)y with a monotonically growing length scale 



after the peak of the dynamic susceptibihty, contrary to reported previously [T^ IT^ l26l |. but in agreement with 
|lll |. We emphasize once more that this result physically makes sense. At times much larger than t* , only very 
rare but very large dynamical domains contribute to the dynamic structure factor, so that spatial correlations are 
weak, but extremely long-ranged. The existence of an ever growing length scale is supported by any model with 
an hydrodynamical limit (such as the phonon or defect models studied here) and is in a sense trivial. The really 
interesting piece of information is the value of this length scale for t = Tq, i.e. when the relevant relaxation processes 
take place. 

We conclude that our numerical data are not inconsistent with a monotonically growing length scale even for 
t > T, although addressing more quantitative issues such as functional form at the growth law and its temperature 
dependence would require quite an important, but certainly worthwhile, numerical effort. 



VIII. CONCLUSION AND FINAL COMMENTS 



Let us summarize the results and the various points made in this rather dense paper. First, we have computed nu- 
merically and analytically, exactly or approximatively, the four-point correlation function designed to characterize non 
trivial cooperative dynamics in glassy systems within several theoretical models: mode-coupling theory, collectively 
rearranging regions, diffusing defects, kinetically constrained models, elastic/plastic deformations. The conclusion is 
that the behaviour of Xii^) is rather rich, with different regimes summarized in the Introduction and in Fig. 1. We 
have computed the early time exponent fi and the peak exponent A for quite a few different models of glass forming 
liquids, and shown that the values of these exponents resulting from these models are quite different, suggesting 
that the detailed study of X4{t,T) should allow one to eliminate or confirm some of the theoretical models for glass 
formation. 

In this spirit, we first simulated some non-cooperative KCMs as the one-spin facilitated FA model in d = 1 and 
d = 3 and the East model. The assumption of point-like defects that diffuse, possibly with an anomalous diffusion 
exponent, gives a good account of the shape of the four-point correlation function and of the four-point susceptibility 
which are in quantitative agreement with the above results for the independent defects model. For strong glasses 
such as Si02, where the relaxation is due to defect diffusion, our results should be quantitative. It would be very 
interesting to reconsider numerical simulations of the dynamics of Si02 under the light of the present paper to check 
in more details that the defect picture is indeed correct in this case (note that our results should enable one to extract, 
in principle, the properties, density and relaxation times of defects from the four-point correlation function). For the 
d = 3 one-spin facilitated FA model, we see clear indications of the interactions between defects as time increases. 
This leads to small deviations of the numerically obtained exponents with respect to those predicted by our analysis 
of the independent defect model, which does not account for interactions between defects. As far as the identification 
of a growing length scale f (<) from numerical data, we have seen that even within this simplified lattice models, this 
can be a rather difficult task. Our results points toward a dynamical correlation length that grows forever and a 
behaviour of S4{k,t) different from the Ornstein-Zernike form but with similar asymptotic behaviour. We leave the 
study of cooperative KCM, for which a more complicated behaviour should occur, for future work. In particular, the 
detailed form of S4{k,t) should contain information about the inner structure of the corresponding defects. 
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We have also analyzed the four-point susceptibility of both a Lennard-Jones system and a soft-sphere system, and 
shown that the initial exponent of the four-point susceptibility is decreasing with the temperature and rather small, 
/i < 1. We have found, perhaps unexpectedly, a reasonable agreement for ^ and A with the predictions of MCT 
but not with other theoretical scenarii, such as simple diffusing defects, strong KCMs or CRR (although this might 
be a question of temperature and time scales, since both CRR and cooperative KCMs are supposed to apply closer 
to the glass transition temperature). Finally we confirm that the extraction of the growth law of ^(t) at a given 
temperature is difficult, and we can only say at this stage that the data is not incompatible with the idea that ^(t) 
grows monotonically, even beyond t = Tq, in the Lennard- Jones system. 

As for further work and perspectives, we think that the following points would be worth investigating. First, it would 
be very interesting to develop a detailed theory of the crossover between the elastic regime described in Section UTTI 
and the Mode-Coupling (3 relaxation regime. Is it possible, in particular, to describe approximately the 'melting' 
of the glass as one approaches the Mode-Coupling transition temperature from below? Second, we only considered 
systems in equilibrium. One in fact expects that the four-point susceptibility also contains very useful information in 
the aging regime (see |13,|H3|)- Detailed predictions in this regime may enable one to probe the mechanisms for slow 
dynamics and the issue of the cooperative length at low temperature in the aging regime |57l |. In particular, the elastic 
contribution should not age whereas the CRR contribution (characterized by the same exponent ^) should exhibit 
some aging, possibly allowing one to separate the two effects. Third, the quantitative study of four-point functions in 
cooperative KCMs where defects have a complex inner structure would be very interesting, since it is clear from the 
present paper that simpler KCMs seem to fail at describing quantitatively X4(^) i^i fragile systems. Fourth, it would 
be interesting to define more complicated correlation functions, for example, a fully general four point function, or 
higher order correlation functions, in order to test in a more stringent way the idea of cooperativity in glassy systems, 
and distinguish systems where the growth of Xi{^) is trivial, such as elastic solids, from those in which a truly non 
trivial cooperativity governs the dynamics. Finally, it seems clear that this issue of cooperativity and its associated 
length scale can only be convincingly settled if long time scales, low temperature regimes can be probed quantitatively 
in experimental systems. We hope that the present paper will motivate ways to directly access four-point functions 
experimentally in glassy systems (see 113); natural candidates for this are colloids Q and granular materials [s^lssj . 
although there might be ways to investigate this question in molecular glasses and spin-glasses as well [s^. 



Fig.|Slwas obtained from the unpublished data of D.R. Rcichman and R.A. Denny. We are very grateful to them 
for providing these results. We thank E. Bertin, O. Dauchot, J. P. Garrahan and D.R. Reichman for discussions. 
G. B. is partially supported by the European Community's Human Potential Programme contracts IIPRN-CT-2002- 
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HPRN-CT-2002-00319 (STIPCO). 
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Appendix A: Dynamics of elastic networks 



A.l The four-point correlation function - Over-damped case 



We will define 6*4 (r,i) for the elastic model defined in the text as: 

Gi{r,t) = (cos g [(?!)(f, i) -^(r,0)] cos (7[(/)(f=0,t) -(?!)(' 
which is equivalent to: 



r = 0,0)]) 



(41) 



Gi{r,t) = -(cosg[</)(r,t)-0(r,O) + 0(r = O,i)-0C 
+ i (cos q [(/)(f, t) - cj)[r, 0) - 0(f = 0, i) + 0(^ 



r = 0,0)] 



f=0,0)] 



Using the fact that the field is Gaussian, we finally find: 




(42) 



where: 



i?(r,t) 



((0(r',t)-0(r,O))(0(r = O,t) 



(/)(f=0,0))) 




(43) 
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Rir,t)^-inty~^/'Fi^) (44) 



Hence, 



with: 



We thus see nTiinediately that G4(r, t) wiU be governed by a 'diffusive' correlation length ^(t) ~ v k^, as expected from 
the structure of the Langevin equation that describes relaxational dynamics. Note that for under-damped dynamics, 
sound waves would change this scaling. 

It is useful to consider the following quantity: 

dl{z) _ f d'^w 



F{z) ^ z^-'[I{^) I{z)]; I{z) = / ^A^e-'^--^ . (45) 



In d = 3, after integrating over dwi, one has 



and: 



^W = i^y^ ^'^du (48) 

Therefore, for z <^ 1, one finds F{z) ~ (47rz)~^ and Rij.t) ~ T/(47rKr), whereas for z ^ 1, F{z) ~ 
(27r'^/^)~^ exp(— z^/4)/z^. Thus, for r <^ ^{t) and nA^t ^ 1, the four-point correlation function behaves as: 

G4(r,t) = /,^(cosh(gl)-l). (49) 
A. 2 The four-point correlation function - Under-damped case 

We have now: 

m^^^ = KA<l>ir,t) (50) 

which has for solutions in the Fourier space: 

0fe(t) =exp(im)(^fc(0), (51) 

with V = {Kjmf-I'^ . We now have: 

([(<^(f, t) - ,^(f, 0)]2) = — y ' = _ y (1 _ cos[Vfct])dfc. (52) 

In (i = 3, we find obviously the same result for jq and G4 as above, but -R(r, t) is now equal to: 

^^^^' = ^ / J^''^''"^' - 

that we write: 

R{f,t)^-[I{r,0)~I{r,t)] (54) 

where: 
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By introducing z = Vtjr and changing variable q = rfc, and also u = cos0 and integrating over u one finds: 

/(r, t)^—\ dqq-\sm{q{\ + z)) + sin(q(l - z))] (56) 

Consider the first term: 

27r /■ 

/(f, t)^—J dqq-^ sin(q(l + z)) (57) 

Changing variable v = q{l + z) directly shows that this integral do not depend on z, as long as (1 + z) is positive. 
This is true for the other integral, which does not depend on z as long as 1 — z is positive. If 1 — z is negative then the 
integral changes sign. Therefore we have that /(r, t) = I{r, 0) if z < 1, and I{r, <) = if z > 1. Therefore R{r, t) = 
if z < 1 and i?(r, t) = ^J^^ when z > 1. The result is very intuitive: when z < 1 the information does not have time 
to travel the distance r and there are no correlation. For z > 1 the two regions are "connected" and one finds the 
free field correlations. Brownian and Ncwtownian dynamics furnish the same correlation for a given r when the time 
diverges, as we expect. Finally, it is straightforward to obtain the result quoted in the text for Xii^)- 



A. 3 The low dimensional case 



We give here, without much details, the results for elastic networks in d = 1 and d = 2. In d = 1, as is well 
known, each particle wanders arbitrary far from its initial position but in an anomalous, sub-diffusing way, as t^/^. 
Correspondingly, the dynamical structure factor decays as a stretched exponential: 

\nCiq,t) ^-qh^/'^ (58) 

K 

Note that the t^^'^ comes from a collective displacement of the cages, and is similar to the anomalous diffusion observed 
for hard spheres in one dimensions, since the latter problem can be mapped onto the Edwards- Wilkinson problem i 
one dimension |58|,i6fl]. We expect that the results obtained here for G4 should also hold for this case as well. In fact, 
this model was recently discussed in the context of a simple d = 1 granular compaction model, see [s^ . 

In d = 2, the displacement grows logarithmically with time, leading to a power-law decay of the dynamical structure 
factor with a q dependent exponent: 

Ciq,t)^t-y y=P^- (59) 

Turning now to X4(0; find that after a short transient, X4(0 grows as t^^'^ in d = 1 and behaves as in d = 2. 



Appendix B: Calculations for the defect model 

In section lvTl we have reduced the computation of Gi(r, t) and Xi{t) to probability distributions of a single random 
walk. In the following we shall show how these quantities can be computed in any spatial dimension. 

Let us call F,^{u) be the probability that a random walk starting in z reaches x for the first time at time u. P^it), 
the probability that a vacancy starts in z at time zero and reaches for the first time x at a time less than reads: 



P^it)^ f F^iu)du (60) 
Jo 

Therefore, we need to calculate F^{u). The trick to do that is writing a linear equation relating F^, that we want 
to compute, to P^{x,t), the probability that a random walk with self diffusion coefficient D, starting in z, is in x at 
time t, which is well known. This linear equation is: 

P'{x,t)=S^Jt.o+ f' F^{u)P''{x,t~u)du . (61) 
Jo 
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By taking Laplace transform (from now on s is the variable conjugated to t and C indicates the Laplace transform) 
we obtain: 

F:it) = c-' (62) 



and 



V £P^(x,s) 



P^it) = f ( ^^T'± J^-' ] {t')dt' (63) 







s CP-\x,s) 



A similar strategy can be used to calculate P^y{t). Indeed the following equality hold 



P^yit) = / Fl^{t')P^{t-t')dt' 

Jo 

Py,^{t) = f Fly{t')PUt-t')dt' 

Jo 



(65) 



where F-^y[t) is the probability that a random walk starting in z at time zero reaches y for the first time at t but 
never touches x at s <t. Therefore, in order to calculate P^^yit) + Py^xit) we need to calculate F^-^{t) + Fy-^[t). It 
is immediate to check that the following equations hold for any choice oi x,z,y 

{ F^{t) = Sx.Jt.o + /o dsFly{s)Fy{t -s) + F^yit) 
\ F; {€) = Sy,Jt,o + Jo dsFl^{s)Fy- (t - s) + Fly (t) 

which implies, again by Laplace transform (z is always different from x and y in the following so we will skip the 
Kronecker deltas), 

^..w+^.w=^- g:i;;:g;;; (67) 

Using the expression for the Fy{s) we get 

F,iy,t)+Fy{x,t)-C -^^j-^—^—^^^-^ (68) 
Furthermore P§{t) = 1 — PJ'(t). Hence we obtain 

£P|W = i-£p.w4(I-|£!||il) ,69) 
Finally, we obtain the expression for 

J-[P.,y[^) + ^y,A^)) CP-{y,s)+CP-{x,s)s[ CP-{x,s)J ^ ' 

An useful way to rewrite this expression is obtained by summing and subtracting the Laplace transform of P^ (t) + 

Pyity- 
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Finally putting together all the different terms we have: 



{n^{t)ny{t)) = exp {~2p„ - 2p^N{t) + 2p^py{t) + p^Git.x - y)) 



(72) 



where Nit) = X^z^a; the average number of distinct sites (minus 1) visited by a random walk during the 

interval of time t and 



G{t,x-y) = C~' 



E 



CP'' (x, s) + CP^ iy,s)l Cpy (x, s) 



s) + CP='{x, s) s CP'^ix, s) 



Since 



{n,{t)f = exp (-2p„ - 2p,N{t))) 

the expression of G4 is 

Giix -y,t) = exp {-2p,, - 2p^N{t)) [exp (2p^P^^(t) + p,,G(<, x - y)) - 1] 



(73) 



(74) 



(75) 



In the following we shall analyze separately the one dimensional case, the three or higher dimensional case and the 
two dimensional case. 



B.l One Dimension 

Consider a symmetric random walk on a one dimensional lattice with lattice spacing a, by Laplace transforming 
the master equation 

dP~(x,t) PHx + a,t) + PHx-a,t)-2PHx,t) , , 

-^r- = 2 

one immediately obtains 

CP%x,s)^ ^^TTY— (77) 

J-^/a 27rC(fc) +s 

where (^{k) = (1 — cosfc). In the continuum limit a ^ 0, (x — y) cx a Dt/2 cx a^, the above integral can be solved 
with the well known result 

CP%x,.s) = -l=e^'^ (78) 

which correspond to the solution of the diffusion equation for a one dimensional Brownian motion with diffusion 
coefficient D, i.e. 



dP d^P , , 



Let us now compute all the functions needed to get 
First 



z^x z^x ^ \ 1 J/ 



where we used equation When t >> 1 we get 
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N{t) = A 



Dt 



Second, using the expression (|73|l of G in terms of CP^{x, s) we get: 



- y^l^r-y 

VD e 



„3/2 -v^l^-B 

* e +1 



Changing variable in the Inverse Laplace transform we get: 



Git) = 4/D</ 



\x-y\ 

V2Dt 



where / ( ) equals 



\x-y\ 
V2Di 



+ioo— 7 



- V27|: 



e -yot 



Finally P^{t) can be computed easily but it is always much smaller than the other terms in the exponential so we 
are going to neglect it. The resulting expression for G4 is: 



6*4(2; ~ y,t) = exp ( —2pi, 



cxp paVDtf 



\x - y\ 
\/2Dt. 



(80) 



Note that the typical time-scale is r = -^^j^ and since we focus on p^, ^ we can rewrite the above expression as: 



G^ix - y,t) ^ cxp 
Integrating over a; — y to get the X4 we find 



V8 



exp 2y^f py 



\x - y\ 



(81) 



Xi{t) = — exp 

Pv 

In particular when t/r <^ 1 we have 



r+00 

^ ) \j2tlT J dx 



exp 



(2^f{x) 



(82) 



Xiit) cx —it/r) 

Pv 



(83) 



The interpretation of this result is that at short times the defects do not intersect and the Xi is just the square of 
the number of average sites visited by a random walk until time t. We'll see that this interpretation is indeed correct 
in any dimension. 

Finally, after some algebra it is possible to obtain from H83|) that Xiit) — ^ exp (^—^y/t/rj at very large times 

(c is a numerical constant). Thus, as found in simulations, the normalized Xi does not go to zero as it happens in 
three dimensions. 



24 



B.2 Three dimension and higher 



Consider a symmetric random walk on a cubic lattice, the general expression for P^(x, s) is 

r jdL ik{x-z) 

where BZ means Brillouin zone and = X)i=i(l ~ cos ki) for a hyper-cubic lattice {ki is the component of k in 
the direction i). Also in this case we consider the continuum limit {x — y)/ ^ Dt/2 oc 0(1) and look for times t much 
larger than one. 

Let us again compute all the needed quantities. First. N{t). In this case for t » 1 we find that 



BZ 



7rC(fc) ) ^ 



Hence N{t) = citD where = /g^ 



Again, we neglect the P^it) term and we focus on G in the continuum limit, for t ^ a we get: 

„ „ _ 1 IbZ (27r)<' Dk-'+s 



\JBZ JTrF)^ Dk 



Changing variable in the Inverse Laplace Transform we get: 



n+ioo-7 I- idi. ik{x-z) px:nffs1 



200 — 7 



BZ (Stt)'' Dk^ + s C^s"^ 



Since we know the Inverse Laplace Transform of the function resulting from the integral over k (it is simply P^{x, t)) 
and each 1/s adds an integral we finally get 



2 I , I , e 2Dti. 



G{t) = C2{Dty du dv 



lo Jo {2DtvY/^ 
where C2 is a numerical constant of order unity. From this expression, we finally obtain 



G4(x - y, <) = exp(-2pi, - 2p„ciDt) 
and the results quoted in the main text. 



exp p,„(c,i5t) duj^ 



(85) 



B.3 Two dimension 

In two dimension things arc bit tricky because of log corrections. Briefly, we obtain that 

G.(.-„,,)^„p(-2-SL^)-c5«/.)'y^ I d.!^ (86) 



with C3 and C4 constants of order unity. Hence, integrating over x — y, we get: 

C3i \ 1 2/,/ n2 1 



In particular when t/r < < 1 we have 

X4(t)(x-f-^)' (88) 

Again, since the number of sites visited in average by a RW in 2D goes like t/\nt, at short times Xi is the square of 
the number of average sites visited until time t. 
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B.4 Density— density correlations 



We now sketch the calculation for the density four point correlation, defined as: 

Gi{x ~y,t) = {{vAt)vAO) - P^)iVy{t)Vy{0) - p')) - {Vx{t)ilM)l 
with i}x{t) = 0, 1 if the site x is empty or occupied at time t, respectively. We start from: 



p 



Using that Ylz P^{x, t) ~ 1 we get 



In the limit p„ ^ we have: 



{Vx{t)Vxml - cxp(-4p,„) (exp(p,P-(x,t)) - 1) 

{vAt)vAO))l = iPvP''ix,t))f 



Similarly we find that 



z,z^x.,y 

= exp (-4p, + 2p,P^x, t) + 2p,py{x, t)) 



Collecting all the pieces together we finally get at leading order in py 

Gi{x-v,t)^2p,py{x,t) 



(89) 



(90) 
(91) 



(92) 



for X ^ y. The interpretation of this eq is that the dynamical correlation between x and y is due to the fact that the 
same vacancy was in x at time and t at time t or vice-versa. Integrating over x — y one finds that at long times 
Xi{t) l/i**/^, showing no interesting structure. 
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